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The transport coefficients of the Anderson model are calculated by ex- 
tending Wilson's numerical renormalization group method to finite temper- 
ature Green's functions. Accurate results for the frequency and tempera- 



^ , ture dependence of the single-particle spectral densities and transport time 

(N , 

2 . r(w,T) are obtained and used to extract the temperature dependence of the 



transport coefficients in the strong correlation limit of the Anderson model. 



0> , Results are obtained for values of the local level position ranging from the 

^ Kondo regime to the mixed valent and empty orbital regimes. The low tem- 

, perature anomalies in the resistivity, p{T), thermopower, S{T), thermal con- 

C , 

O . ductivity k{T) and Hall coefficient, Rh{T), are discussed in terms of the 

CJ 

behaviour of the spectral densities. All quantities exhibit the expected Fermi 
liquid behaviour at low temperature, p{T) = pq{1 — c{T/Tk)'^), S{T) ~ 7T, 
K{T)/a T = 1 + (i{T/TK?, Rh{T) = -Roo{l - S{T/TKf). Analytic results 
based on Fermi liquid theory are derived here for the first time for /? and the 
numerical results are shown to be consistent with this coefficient. The range of 
temperatures over which universal behaviour extends is also discussed. Scat- 
tering of conduction electrons in higher, / > 0, angular momentum channels is 
also considered and an expression is derived for the corresponding transport 
time and used to discuss the infiuence of the interference terms between the 
resonant / = and non-resonant / = 1 channels on the transport properties. 
The presence of non-resonant scattering is shown to be particularly impor- 
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tant for the thermopower at half-filling, where the sign of the thermopower 
can depend sensitively on the non-resonant phase shift. Finally the relation 
of the results to experiment is discussed. 

INTRODUCTION 

In this paper we present accurate results for the transport coefficients of the Ander- 
son model obtained by extending the numerical renormalization group method to the 
calculation of finite temperature Green's functions. The Anderson model has been used 
extensively to interpret the properties of dilute magnetic alloys and a number of the local 
properties of heavy fermion compounds. It is also of use in discussing the properties of 
concentrated Kondo systems in cases where alloying or disorder inhibit the coherence ef- 
fects between the magnetic ions. Although many of the properties of this model are now 
well understood 0, the temperature dependence of the transport coefficients has proved 
to be particularly difficult to calculate reliably. The transport coefficients require accurate 
expressions for both the temperature and frequency dependence of the impurity Green's 
function, a quantity which is difficult to calculate in the strong correlation limit of the 
Anderson model. Here we calculate this quantity from Wilson's numerical renormalization 
group method |l|] which is non-perturbative and is therefore accurate for arbitrarily large 
Coulomb interactions. This method has played a crucial role in forming our current un- 
derstanding of the Anderson model. It was first applied to the Kondo model by Wilson 
HI] and subsequently to the Anderson model by Krishnamurthy et al 0. The two models 
are related via the Schrieffer- Wolff transformation [Q, with the Kondo model describing the 
low-energy physics of the Anderson model in the strong correlation Kondo regime. The 
application of the numerical renormalization group to these models yielded the elementary 
excitations, thermodynamics, fixed points and effective Hamiltonians around the fixed points 
Dynamic and transport properties were not calculated. The calculations showed that 
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the Kondo model has two fixed points which characterize its physical properties. The local 
moment fixed point which describes the high temperature regime and in which the con- 
duction electrons couple weakly to the impurity moment and a strong coupling fixed point 
which describes the low temperature regime in which the impurity moment is quenched and 
the excitations are those of a local Fermi liquid. The parameter space of the Anderson 
model is larger and the calculations showed that in addition to the local moment and strong 
coupling fixed points there are two additional fixed points. The most important of these 
is the valence fiuctuation fixed point and is characteristic of the asymmetric model. In the 
valence fiuctuation regime, charge fiuctuations become important and the properties corre- 
spond to a model with a strongly renormalized temperature dependent resonant level 0. 
In this paper we describe the transport coefficients and their relation to the corresponding 
spectral densities in these different regimes, pointing out the characteristic features which 
arise in each case. 

In contrast to thermodynamic properties, which have been obtained exactly by the Bethe 
Ansatz |^,^ and numerical renormalization group the calculation of transport and 

dynamic properties have relied on approximate methods. Finite order perturbation in U 
calculations give accurate results in the Fermi liquid regime for spectral densities and ther- 
modynamic properties up to U/tiA 2.5, where U is the local Coulomb repulsion and A 
is the resonant level width in the Anderson model 0,^]. However, in the local moment 
regime U >> A for T » where properties depend on logarithmic terms this approach 
breaks down. Quantum Monte Carlo approaches |TO| , p!I| become increasingly less accurate 
for larger values of U and lower temperatures. So far transport coefficients via this method 



have been calculated for just the symmetric Anderson model and for U/irA < 3 |T2[. Trans- 
port coefficients have also been obtained for larger degeneracies of the local level via the 
non-crossing approximation |1^,|T^. A problem with this method is that it fails to satisfy 



the Fermi liquid relations at zero temperature [|T5|. In a finite magnetic field the transport 
coefficients for the Kondo model were discussed in []T^ on the basis of the Nagaoka inte- 
gral equations. Recently a comprehensive and highly accurate approach to the calculation 



of dynamic properties of magnetic impurity models has been developed by extending the 
numerical renormalization group approach [p!7|-p3|. This overcomes the above mentioned dif- 
ficulties with the approximate methods. Accurate results in all regimes have been obtained 
for single-particle spectral densities at both zero [p!^-p2| and finite temperature p3| , pl 



These satisfy all the sum rules and Fermi liquid relations. In the next section we introduce 
the Anderson model including terms which model the scattering of conduction electrons in 
higher (/ > 0) angular momentum channels. The transport coefficients are defined in terms 
of the transport time for conduction electrons scattering in both the the resonant and non- 
resonant channels (the transport time incorporating non-resonant scattering of conduction 
electrons is derived in Appendix 1). The numerical renormalization group and its use in 
extracting finite temperature Greens's functions and spectral densities is then described. We 
also give analytic calculations for the low temperature behaviour of the transport coefficients 
based on Fermi liquid theory. These are used as a check on the accuracy of the numeri- 
cal results. Finally we present the conclusions and indicate the relevance of the results to 
experiment. 



MODEL, TRANSPORT PROPERTIES AND METHOD 

The model 



The Anderson model including non-resonant scattering of conduction electrons in / > 
channels is given by the following Hamiltonian 

H = Himp + Hhyb + He, (1) 

+/ 

Himp = eocLcoa + Uno^noi + (^l^m'^lm, 

<y l>0 m=-l 

+1 

Hhyb = Y^{VQkc{^CQ^ + H.C.) + Y i^^^kcl^^lm + H.C.), 

kcr ka l>0 m=—l 

He = Y^kci^^ka- 

ka 
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The first term Himp represents the impurity and is coupled to the conduction electrons He 
via the hybridization term Hhyh. The scattering of conduction electrons in Z > channels is 
modelled by including uncorrelated levels, e/, hybridizing with the conduction electrons. This 
is equivalent to taking into account phase shifts rji with / > for the conduction electrons 
in addition to the usual / = resonant phase shift. The charge neutrality condition requires 
that these phase shifts satisfy the Friedel sum rule [§] 

z = j2m+i)-, (2) 

1=0 ^ 
where Z is the excess charge on the impurity. 

The many body effects arise from the strong Coulomb repulsion between the electrons 
in the impurity / = level. 



Transport coefficients 

Assuming that the conduction electrons scatter incoherently from a small concentration, 
rii << 1, of magnetic impurities, linear response theory allows the transport coefficients to 



be expressed in terms of the transport integrals ||28 

'•+°° / df{uj ' 



-'ml 



T {uj){uj — n)"^duj, (3) 



where fi is the chemical potential, and t{u!,T) is the transport time. The resistivity, p(T), 
thermoelectric power, S{T), thermal conductivity k(T) and Hall coefficient Rh{T) are given 
in terms of these by 



1 Ln 


1 /^r(^,T)(-f 


) duj 


e 




e 


T It{uj,T) 







k{T) = Ul,,-^], (6) 



T L, 
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Rh(T) — —RoojY', 



(7) 



where = ni\e\c. In the absence of non-resonant scattering, the transport time, tq, (see 
Appendix 1 for constant factors such as rii) is given by 
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Apo(c^,T), 



(8) 



To(a;,T) 

where A is the resonant level width and po is the single-particle spectral density. The latter 
is given in terms of the resonant level Green's function Go{uj,T) =« Co„;cl^ » and 
self-energy T,{uj,T) = T,^ + iT,^ by 



Po(w,2") 



(A - E^cv)) 



(9) 



The transport time in the presence of non-resonant scattering is derived from the Kubo 
formula for the conductivity in Appendix 1. In contrast to the case of resonant scattering 
only (see e.g. [0), the vertex corrections for the current-current correlation functions enter- 
ing the expressions for the transport coefficients are finite when non-resonant scattering is 
included. The resulting expression for the transport time after inclusion of vertex corrections 
is 
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T(a;,T) ro(u;,T) 



Re Go{u,T) . 
cos 2r]i - - — — — — sm 2r]i 
Im Gq[uj, I ) 



+ Pn 



(10) 



The effects of non-resonant scattering are primarily contained in the factor 



o Re Go(u,T) ■ o 

cos2r7i - ^^^»;(^y^ sm2?/ij 
nels. The non-resonant, / 7^ 0, phase-shifts in ( [lOl) are taken to be constants, defined by 
the screening charge in respective channels. Thus, 



which is due to interference between the / = and / = 1 chan- 



Pn 



mkf 



sin^ r]i + J2^ sin2(?7, - ?7,_i) 
i>i 



(11) 



In the T = limit (but with / 7^ 0) the transport coefficients calculated with (^) - ([TlD 
reduce to standard phase-shift expressions, while in the limit rji^o — > we recover the 
usual many-body expression for the single-channel transport time (S). 
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From the expressions for t{uj,T) and to{uj, T) we see that in order to evaluate the trans- 
port coefficients we require an accurate expression for the frequency and temperature depen- 
dent resonant Green's function Go{uj, T). We obtain this from the numerical renormalization 
group method as described below. 

The numerical renormalization group method 

The numerical renormalization group for the Kondo and Anderson impurity models is 
described in , where it was used to obtain the thermodynamic properties. Here we give a 
description of the method and its use in calculating finite temperature Green's functions and 
specifically the local Green's function Go{uj,T) required for the transport time. The central 
idea in the numerical renormalization group is the importance of including all energy or 
length scales. The Hamiltonian (|I]) contains conduction electron states of all energies from 
the band edge D down to zero energy and states from each energy scale contribute to the 
impurity properties. To take into account these states Wilson introduced a logarithmic 
discretization of the conduction band about the Fermi level so that all energy scales were 
represented, with the greatest resolution at low energies where the many body effects are 
most important. As shown in Appendix 2 this logarithmic discretization approximation 
results in the following discrete Anderson model for the resonant channel, 

H = UrriN^^^il + A-i)A-(^-i)/2^^^ 

N-l 

= A(^-'y^[Ho + Hny, + E ^""'^Ufl+iJu. + fLUi.)i (12) 

n=0 

where Hq = eocl^Coa- + UnQi^noi is the resonant part of the impurity, Hhyt = Vo{fQ„Co^ + H.c.) 
couples the impurity to a local conduction electron orbital /oo-|0 >, and the last term de- 
scribes the remaining conduction electron orbitals whose wavefunctions have a large overlap 
with the impurity. The conduction electron orbitals neglected in the above discrete approxi- 
mation to the full Anderson model have their wavefunctions localized away from the impurity 
site and have negligible contributions to the impurity properties (see Appendix 1 and for 
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further details). The parameter A > 1 describes the discretization of the conduction band. 
The above discretized Hamiltonian is in the form of a semi-infinite hnear chain and can be 
iteratively diagonahzed for increasing chain lengths to obtain the eigenvalues, , and 

N — 1 

eigenvectors, \p >n , on successively lower energy scales un ~ -DA where un is the 
lowest scale of (see Appendix 2 for details). From the eigenvalues, the thermodynamic 
properties are obtained at a corresponding sequence of temperatures T/v ~ cutv/^b- The 
details are given in and in Appendix 2. Here we show how the local Green's functions 
can be extracted on successively lower energy scales. 

Consider the Green's function Gq{ijJ,T) =<< cqo-; cj^ >>• Using Hn, the N'th shell 
Green's function Gq {uj,T) and associated spectral density Pq {uj,T) are given by 

-l'.,.,,. (13) 
(a,, T) = -\- J: |M^,f ( e-^«? + e-^^? )S{^ - E« + (14) 

NyH ) pp' 

Here Z^iP) is the partition function for the N'th cluster, and M^, =< p\coa\p' > are the 
many-body matrix elements of the local operator cqo-- The latter can be evaluated recursively 
in a similar way to the evaluation of the matrix elements < p\fNa\(l > in ( |B14| ). Using the 
unitary transformation ( [B13|) we obtain 

M^p, = < p\coa\p' >N= U*NiP^ ri)UN{q, r'i') < i\ < r\co^\r' > \i' > 

r,i r' ,i' 

= Y.llU*N(P^rt)Ur,{q.r't')5,,,M^^;;\ (15) 

Hence the matrix elements M^pi can be evaluated recursively from a knowledge of the eigen- 
states of the N'th cluster, Un{p, ri), and the matrix elements, M^r^, of the previous cluster. 
For a given cluster size A^ the Hamiltonian only describes excitations in a limited range 
of width Kujn above the lowest scale un present in H^, due to the truncation of the spec- 
trum as described in Appendix 1. At T = the spectral density is evaluated at u ^ 2un- 
Calculating the spectral density at energies much smaller than this using Hn is not justified, 
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since information on lower energies is obtained in subsequent iterations, whereas calculat- 
ing the spectral densities at much higher energies than this might introduce errors due to 



the truncation of the spectrum on the high energy side. In evaluating (|T4|) the delta func- 
tions are broadened with Gaussians of width a at of order ojn appropriate to the cluster 
size. The broadening parameter a^r used within each energy shell is varied continuously so 
that there is no discontinuity in going between successive iterations. The small remaining 
asymmetry in the spectral features due to the larger broadening parameter at the higher ex- 
citations should vanish in the limit A ^ 1. The procedure for calculating finite temperature 
Green's functions is slightly more complicated. The shell Green's function Gq{uj,T) and 
spectral density Pq{uj,T) are only reliable for excitations u ^ 2uji\f and for temperatures 
/c^T << 2co'7v. For temperatures ksT » 2un excited states not contained in Hj^ would 
be important, whilst for temperatures ksT of order 2ujn- transitions between excited states 
would need to be known accurately. In principle these are known from subsequent itera- 
tions, but are not contained with sufficient accuracy in H^. The only transitions known with 
sufficient accuracy in Hj^r are the groundstate to excited state transitions with excitation 
energies ~ 2co'Ar which is the natural energy scale of this cluste-r^. As long as ksT « 2un it 
is not necessary to know the lower excitations since their contribution to the Green's func- 
tion and spectral density for the energies u = 2ujn under consideration will be negligible 
(note the delta function in (0)). From this discussion it follows that the spectral density 
for temperature T can be calculated from the shell spectral densities Pq{uj, T) at frequencies 
uj ^ 2uJi for i = 1, 2, . . . , M until 2ujm becomes of order T. To calculate the spectral density 
at temperature T and for frequencies below 2ujm, a smaller cluster is used. This is done 
because when T is larger than the frequency at which the spectral density is being evaluated, 
it is the excited states of order T contained in previous clusters which are important and 
not the excitations very much below T. The procedure outlined here requires storing all the 



^The term cluster is suggested by the notation although the calculations are in k-space 
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matrix elements for the Green's functions for each cluster size, since smaller clusters may 
be required in subsequent iterations. 

In the absence of non-resonant scattering we can calculate the transport time tq{uj,T) 
directly from the single-particle spectral density po{uj,T) using (||) and hence the transport 



real part of Go{u!,T) and use (|T0|) and (§-0) to calculate the transport coefficients. In the 
next section we present the numerical results obtained using this procedure. 



Before presenting numerical results we outline some analytic results on the low temper- 
ature behaviour of the transport coefficients obtained by using Fermi liquid theory. These 
serve as an independent check on the accuracy of the numerical renormalization group 
method. In the following analytic calculations we restrict ourselves to the case of resonant 
scattering only so the transport time t{uj,T) = To{uj,T). 

To extract the low temperature behaviour of the transport coefficients we use the Som- 
merfeld expansion. In the transport integrals @ the factor ) for temperature T pro- 
vides an energy cut off outside the Fermi window Ic^l < T. The functional form of Gq{uj,T) 
also has an energy scale << ksT^ over which t{uj,T) is a a slowly varying function, 
where Tq = Tk for the Kondo regime and ksTo = A, for the mixed valent regime. As a 
consequence, for temperatures T « Tq, we can apply the Sommerfeld expansion. 



The quantity h{u!,T) entering the above expression is t''{u!,T)u"^ for the integral Lmi- For 
the Kondo problem the transport time t[u!,T) is a strong function of temperature at the 



coefficients In the presence of non-resonant scattering, we evaluate in addition the 



Low temperature results 





(16) 
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Fermi level, so in the Sommerfeld expansion in addition to the temperature dependence 
coming from the fermi function there is also the explicit temperature dependence of the 
transport time at the Fermi level. This is taken into account by expanding T{eF,T) in a 
Taylor series 

r(e^,T)=r(e^,0) + ^ + O(T^) (17) 

In the Sommerfeld expansions for the resistivity, thermopower, thermal conductivity 
and Hall coefficient the derivatives which enter are ^f^, and ^^jk- Since to(uj,T) = 

1/ Apo{u!,T) these derivatives are given by 



1 dr 


1 9po 




Po ' 


1 d^T 


^ (po ) 


T duP' 




1 


/ 1 5Vo 


r(9T2 





(18) 
(19) 
(20) 



Using d^) we obtain for the energy derivatives oX uj = ep and T = 0, 



l/5po^ 
Po \ . 



37cot(7™o/2), (21) 



1 9Vo\ 1 , /27rpo(eF,0)^' 



Po / A \ 9a;2 



cot (7rr2o/2) 



o . . .^ / 1 1 /a^S^X COt(7™o/2) /a^S^W 

+ 2.Apo(e., 0) (-^ + ^ j j j , (22) 

where A = zlS. and z = (I — is the wavefunction renormalization constant. In 

deriving the above we have used (|^) and the Fermi liquid properties of the self-energy 
T,^{uj, 0) ~ and 0) ~ (cj — e^)^. The quantity hq is the local level occupancy and 7 

is the linear coefficient of the low temperature specific heat. The latter is given by ||^ 

^ 3 ''e~2 + A2' ^ ' 

where 60 = z(eo + S^(ei7',0)). In the Kondo regime z « 1 and A ~ ksTx « A which 
leads to a large enhancement of 7 since from (^) 7 ~ I/Tk- The Kondo temperature, Tk, 
is given by |^ 



11 



(24) 



For the thermopower we require only the first energy derivative of po to extract the 
leading term. We have from (|^) and (|16|), 



S{T) 



1 



\e\T 



Tier) 



2 1 



(25) 



which on using (|T8|,21) reduces to the Fermi hquid resuh, |33.35 



S{T) = ^ cot(7™o/2) + 0{T^) 



(26) 



and can be checked by evaluating S, 7 and no within the numerical renormalization group 
technique |^ . 



To extract the other transport coefficients to lowest order in T requires the second deriva- 
tive with respect temperature of tq and hence an expression for the self-energy correct to 
order T^. This is difficult to obtain analytically in the general case so we restrict ourselves 
to just the symmetric case. The self-energy for the symmetric model correct to order uj"^ 
and is given by [|T],|2| 



S(cj, T) = S^(cj, T) + tJ:\uj, T) = S^(cj, 0) 







2A(7rA)2 

where Fq is the irreducible vertex function evaluated at zero frequency 
we obtain the temperature derivative, 



(27) 



Using 



9T2 



(28) 



'.=e,,T=o KPodT^J vrA ((6o + S«)2 + A^, 
We now use these results to first check the exact result for the coefficient of the resistivity 
first derived by Nozieres [^] and Yamada [^,0 and then derive the exact coefficient for the 
coefficient of {k,{T) /T) / {k{T) /T)t^o- Restricting the discussion to the Kondo regime 
where the charge susceptibility is zero so Fq = vr A / z, and the particle- hole symmetric model 
where eo + = 0, we find from (E^, p^) 
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and 



TdT^ 



(29) 



Hence from (^, p!6| , p!7D and (|2^) we have for the resistivity, 

1 



piT) 



lrr2 1 aV I n^kBTf 

-1 rm — TTToT^ "I 



1 



6 T(eF,0)a^ 



e2r(ei.,0) 
p(0)(l-c 



2 A2 



1^ 
2 A2 



rp2 + C)(T4 



T 



K 



(30) 



where p(0) = l/eV(eF,0) and c = TT^k%T^/A^ = nyiQ = 6.088 and we have used the 
result for the symmetric case vrA = AkBT^ [^,0. The above result for c is the exact result 
obtained by Nozieres |30]. Without the contribution from (E^) the coefficient c would have 



been wrong by a factor of 2. Away from the symmetric case the result of Nozieres will be 
modified due terms in (|TH|-|T^) which depend on the occupation number no and also due to 
additional temperature dependences coming from the occupation number and the imaginary 
part of the self-energy. These are difficult to estimate analytically, but our numerical results 
to be presented in the next section give some indication of their size. Most terms depend 
on cot(7r?2o/2) and are expected to be small in the Kondo regime no ~ 1. 

We now prove an exact Fermi liquid result for the thermal conductivity. Using 
\e\S(T)/T = 7r7 cot(7rno/2) we have 

<T) 1 



T 



The expansion for L21 is 



2^2 



L21 + T^l cot(7rno/2)L 
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(31) 



+ 



720 



(kBTY 



'd\uj^T{uj)y 



TT^klT^ 



r(ep,0) 



3 _a 



2n^kl 



1 + 



27r2A;| 



+ 



T=0. 



T=0 



60 



{kBT) 



+ 0{T 



6^ 



+ 0(T^ 



20 r I du^ 



0{T' 



■ (32) 
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Substituting this into the expression for the thermal conductivity, and using the above results 
for the derivatives we obtain for the symmetric case in the Kondo limit, 



Ti^k\T{tFi 0)/3. In the asymmetric case there will be additional contributions to the thermal 
conductivity which will modify the coefficient (3. As discussed earlier for the resistivity these 
will be difficult to estimate but are expected to be small in the Kondo regime, so we expect 
the coefficient (3 to remain close to its symmetric value 137r^/80. 

The low temperature behaviour of the Hall coefficient can also be extracted. From the 
Sommerfeld expansion it can be shown that the order terms from the energy derivatives 
cancel and that the term results from the temperature dependence of r at the Fermi 
level, yielding in the symmetric case Rh(T) ~ — i?oo(l — ^(T/Tk)"^) where the coefficient S 
depends on second derivatives and a fourth derivative {d^r / duj'^dT'^)^^=^p^T=o- We do not 
evaluate this since as discussed in the section on numerical results, there is a much larger 
contribution to the Hall effect due to skew scattering of electrons which should be included 
for a discussion of the Hall coefficient of magnetic impurities [^,^]. 

To summarize the analytic calculations, we see that in Fermi liquid theory in the Kondo 
regime that the transport coefficients all scale with at low temperature. Specifically we 
have 




(33) 



where 13 = 26c/10 = 1371^/80 = 15.83 and the linear 




p(T) = p(o)(i-c(T/r^n, 



(34) 



S{T) 



cot{r]o{ep)) ~ (T/Tk), 



(35) 



e 



{n{T)/T)/iniT)/T), = l + P{T/TKr 



(36) 



Rh{T)IR, 



'OO 



{l-5{T/TKf). 



(37) 
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where the constants c, a,/5,7 and S depend on r and its derivatives at the Fermi level. For 
the symmetric model one finds, 

c = — = 6.088, (38) 
a = i^)o = ^r{e„0), (39) 
p = (13/80)7r^ = 15.83, (40) 

Corresponding expression for arbitrary U can be derived from a renormalized perturbation 
theory in terms of A and a renormalized interaction U = z'^Tq p7 |. 

In the asymmetric case for the Kondo regime no ~ 1, there should be negligible cor- 
rections to these coefficients, the asymmetry in this regime should only affect the value of 
Tk. 

The functional form of the low temperature transport coefficients in the Kondo regime 
(P^-|37D and the coefficients of the leading order terms (p8|^TD in the symmetric case provide 
a check on the accuracy of the numerical renormalization group results presented in the next 
section. 



NUMERICAL RESULTS 

The numerical results presented here are for the strong correlation limit of the Anderson 
model, i.e. U/nA » 1. We take U/ttA = 4, unless otherwise indicated. The local 
level position takes values eo/A = —U/2A = —2n corresponding to the symmetric case 
and Co/A = — 4, — 3, — 2, — 1, 0, +1 corresponding to the asjTumetric case. We are mainly 
interested in universal effects independent of the band structure, so all parameters are small 
relative to the conduction electron half-bandwidth D = 1. Specifically we used A = O.OID. 
For eo/A = —4, —3, —2 the parameters describe the Kondo regime, whilst for eo/A = —1, 
and eo/A = +1 they describe mixed valent and empty orbital regimes respectively. In 
discussing the effects of non-resonant scattering on the transport properties we took a small 
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/ = 1 phase shift, r]i = iO.Olvr, ±0.027r, . . . , iO.OSvr. The discretization parameter A can be 
taken as low as 1.5 without encountering large errors due to the truncation of high energy 
states. The use of such a small value of A does introduce errors into the spectral densities 
for extremely low energies. These are noticeable in the symmetric case for example, where 
for sufficiently low energies u « ksTx, the symmetry of the spectral density is violated 
slightly. A larger value of A allowing the lower energies to be reached in fewer iterations could 
be used to avoid this problem, however this was not necessary in the present calculations. 



Thermodynamics 

The thermodynamic properties of the Anderson model have been discussed in where 
the static susceptibility was calculated over the whole temperature range from the band edge 
down to temperatures T « Tx- The specific heat was calculated in the strong coupling 
regime T << by perturbation theory on the effective Hamiltonian around the strong 
coupling fixed point [0. Here we briefiy present the results for the local level occupancy 
uqIT) and specific heat C{T) over the whole temperature regime of interest (calculations of 



C(T) for the screened Anderson model where presented in |21 ). 



Local level occupancy, no{T) 

The local level occupancy rioiT) can be calculated at a sequence of temperatures T/v, = 
1, 2, . . . from the average electron number on the impurity and is shown in Fig. |l]for several 
values of the local level position ranging from the Kondo regime (eo/A = —4, —3, —2) to the 
mixed valent (cq/A = —1,0) and empty orbital regime (cq/A = +1). The only significant 
temperature variation occurs on the valence fiuctuation temperature scale fc^T = A = 
10~^D. On the Kondo temperture scale, (listed in Table the variation of uq is 
insignificant indicating that the Kondo resonance has negligible weight in agreement with 
the spectral densities presented below. The values of nQ{T = 0) are shown in Table |. 
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Specific heat 



The specific heat, shown in Fig. 0, exhibits one or two distinct peaks depending on the 
position of the local level. In the Kondo regime there are two distinct peaks, one at T 
and another at T of the order of A. The low energy peak, which is observed in many 
systems, is associated with the entropy of the magnetic impurity so the area under this 
peak is given by kBln{2) for the present case where the impurity has S = 1/2. The higher 
peak is associated with the charge degrees of freedom. The same behaviour is found in the 
Bethe-Ansatz solution |^|. In Fig. ^we show the ratio C{T)/T. The linear in T behaviour 
characteristic of a Fermi liquid is clear and we tabulate the linear coefficient in Table |. 
Finally, Fig. ^ shows that the specific heat is a universal function of T/Tk in the Kondo 
regime. Universality for this quantity extends up to T ^ 5Tx for the case eg/A = —4. From 
the figure it seems very likely that in the Kondo regime and provided the local level is very 
far from the Fermi level the universal behaviour should extend over the whole range lOT^ 
of the low low energy peak. 

In the mixed valent regime there is a peak slightly below the valence fiuctuation tem- 
perature scale ksT = A and a shoulder or small peak at higher temperatures. The same 
behaviour is found in the empty orbital case. The same behaviour is found also from the 
Bethe-Ansatz solution . 

The value of the linear coefficient of specific heat has been checked in the case of f/ = 
and found to be in good agreement with the exact result. It has also been checked 
independently for both U = and finite U |^8| by using the Fermi liquid relation, 

^""^ +Xc = 7r^7, (42) 



where Xs and Xc are the spin and charge susceptibilities and can be obtained as thermo- 
dynamic quantities or from the corresponding dynamic susceptibilities. This Fermi liquid 
relation has been shown to be satisfied within a few percent and provides an independent 
test on the accuracy of 7 from the present thermodynamic calculation pS] 
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Spectral densities 



The spectral densities in the Kondo regime are shown in Fig. |^-|^ for zero and finite 
temperatures. At low temperatures T « the spectral density is characterized by three 
resonances: the Kondo resonance at the Fermi level of width Tk and negligible weight and 
two atomic-like resonances on either side close to the bare excitation energies eo and eo + U. 
The latter contain most of the spectral weight and have widths of order A. The Kondo 
resonance lies at the Fermi level in the symmetric case ep = —U/2 and above the Fermi level 
in the asymmetric case cq > —U /2 (see inset to Fig. The width of the Kondo resonance is 
approximately Tk, where Tk is defined in (^) and shows the correct exponential dependence 
on eo and U (inset to Fig. ^and [^). On increasing the temperature, the Kondo resonance 
decreases rapidly in height and eventually disappears completely for T >> Tk (Fig. 
The atomic-like peaks remain temperature independent until ksT ^ A. For fc^T > A 
they acquire some temperature dependence, and in the asymmetric case. Fig. 0, there is a 
transfer of spectral weight from the lower to the upper peak with increasing temperature. In 
the symmetric case. Fig. no spectral weight can be transferred but instead the two peaks 
broaden with increasing temperature. The positions of these high energy peaks changes 
little with temperature. Qualitatively similar results for the temperature dependence of the 
spectral densities have been obtained by perturbative methods for fZ/vrA < 2.5 [0]. The 
spectral densities in the mixed valent and empty orbital regimes are shown in Fig. [sHTOl for 
zero and finite temperatures. In these regimes the spectral density at low temperatures, 
ksT << A, is characterized by two resonances (see Fig.pHTOl): one at cq of width approxi- 
mately A and a much broader resonance at eo + U carrying very much less weight. In the 
mixed valent regime, eo/A ^ —1, the bare level is renormalized by the interactions to lie 
above the Fermi level e~o > 0. The bare level is also renormalized to higher energies in the 
empty orbital case , i.e. e~o > eo > A. On increasing temperature, the main effect in the 
mixed valence case is a strong temperature dependence of the resonance at eo on a scale of 
order A (see Fig. p|-p!0|). This resonance broadens and decreases in height with increasing 
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temperature, but does not completely disappear at high temperature. In addition it always 
remains above the Fermi level. The resonance at e~o + U remains distinct from that eo up to 
at least fc^T ^ 6 A. Similar behaviour occurs in the empty orbital regime but the relevant 
energy scale for the temperature dependence in this case is set by eo- 

A good test on the accuracy of the spectral densities is provided by the Friedel sum 
rule which relates the spectral density at the Fermi level to the T = occupation number 
no(r = 0): 

Poiu = eF,T = 0) = ^ sin2(7rno/2). (43) 

Using uqIT = 0) calculated from the partition function gives agreement to within 3% in all 
parameter regimes as shown in Table |. The values of uqIT = 0) calculated from the spectral 
densities are within 3% of those calculated from the partition function, but the former are 
expected to be essentially exact since they are obtained from the low energy part of the 
spectrum {nQ{T = 0) = limj-^ < '^o >Tm) where the method gives the highest resolution. 
It should be noted that the value of uq calculated from the spectral density is largely due to 
the high energy features at eo, where the method has less resolution due to the logarithmic 
discretization of the conduction band. However, the agreement to within 3% between the 
values of no calculated by these two different ways provides additional evidence that method 
can give accurate results for high energy spectral features as well as low energy ones. The 
high energy features have a small asymmetry due the broadening procedure used for the 
delta functions in the spectral density which should vanish in the limit A — > 1 (as described 
earlier) . 



Transport coefficients 

Resistivity 



Results for the resistivity are shown in Fig. 11 for several values of the local level 



position ranging from the Kondo regime (eo/A = —^,—4, —3,— 2) to the mixed valent 
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(eo/A = —1,0) and empty orbital (eo/A = +1) regimes. The behaviour of the resistivity is 
quahtatively similar in all regimes with a monotonia increase with decreasing temperature. 
At T = it reaches its limit p(0) = l/e^ro(ei7', 0). Since tq ~ 1/po "we have from the Friedel 
sum rule for the spectral density that the T = resistivity p(0) ~ sin^(7™o/2). This is 
satisfied for all cases shown in Fig. |TT] with the same accuracy that the Friedel sum rule was 
satisfied for the spectral density (discussed earlier and summarized in Table |l|). 

In the Kondo regime at low temperature T << Tk the resistivity is given by the exact 



Fermi liquid result [l30|-p2 



p(T)=p(0)|l-c(^)'|, (44) 

4 

where c = jg = 6.088 with corrections which increase with increasing asymmetry, as indi- 
cated in the previous section. The rescaled resistivity in the Kondo regime is shown in Fig. 



T3| . Universal behaviour is found in the range < T < 5Tk- The inset shows the expected 
Fermi liquid behaviour. In the Kondo regime the values of the coefficient extracted 
from a least squares fit in the region < T < 0.1 Tr- are within about 8% of the exact result 
c = 6.088 (see Table |). The values of c increase systematically with increasing asymmetry 
and we believe this is in part due to the correction terms discussed previously. However the 
increase in c is small, less than 8% in going from the symmetric to the eo = — 2A case, which 
indicates that the correction terms discussed earlier are also small in the Kondo regime. 
At high temperatures T » T^, the resistivity is well described by the Hamann result 
for the sum over parquet diagrams, 

m = P^(i Ht/Tkh) \ ( . 

^ 2 [ [{HT/Tkh)' + 7^'S{S+1)Y/^)- ^ ^ 



Fig. shows the resistivity for the symmetric case together with a fit of the Hamann result 
with Tkh = Tk/^-2- The fit is very good in the range <T< lOT^ but ( ^5]) fails to give 
the correct Fermi liquid behaviour at low temperature. 



The effect of non-resonant scattering on the resistivity is shown in Fig. 11^ for the sym 



metric case. The resistivity decreases slightly at low temperature for increasing absolute 
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values of the non-resonant phase shift \rii\. The overall behaviour, however, for small val- 
ues of Ir^il is the same as for resonant scattering only. The decrease with increasing |?7i| is 
due to the induced change in the local charge 6n away from uq = 1. Since the resistivity 
attains its maximum for no = 1, any change away from this leads to a smaller resistivity. 
In the asymmetric case, the effect of non-resonant scattering is to decrease the resistivity 
for increasingly negative values of t]i and to increase it for increasingly positive values of t]i. 
This effect is also due to the induced change in the local charge due to a repulsive {rji > 0) 
or attractive (rji < 0) local potential. 

To summarize, the resistivity is described accurately over the whole range of temperatures 
from the logarithmically dominated region T » Tx through the crossover region T ~ 
and into the low temperature Fermi liquid region T « T^- 



Thermopower 



The thermopower is shown in Fig. |T5|-0 for the case of resonant scattering only. 
The effect of non-resonant scattering on the thermopower is discussed below. In Fig. ^ 
we show the low temperature Fermi liquid behaviour of the thermopower given by (p6|), 
S{T)/'~fT = (tt/c) cot(7rno/2), which relates the linear coefficients of the thermopower and 
specific heat to the resonant level occupancy. By extracting 7 from the second derivative 
— limy >o d"^ Fimp{T) / dT"^ of the impurity Free energy and no from the spectral density or 



partition function we have shown elsewhere |24| that this Fermi liquid relation is satisfied 



to a high degree of accuracy in all parameter regime^. 



From Fig. p!q-[T7| it can be seen that the thermopower has different characteristic be- 
haviour in the Kondo and mixed valent regimes. The thermopower for the symmetric case, 
€0 = —U/2, is identically zero due to particle-hole symmetry (see later for a discussion of the 



In the inset to Fig. 4 of |24] the curves S{T)/T were labelled in the reverse order to the correct 
one shown here for different local level positions 
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effect of the non-resonant terms in this case). For the asymmetric case in the Kondo regime 
(eg/A = —4, —3, —2) it exhibits a low temperature maximum at T ^ and then changes 
sign for T > Tk before reaching a broad minimum at high temperature fc^T ^ A (Fig. ^6|). 
In the mixed valent regime (cq/A = —1,0) there is again a maximum, at ksT ^ A/3, 
but there is no sign change at higher temperature and there is only a shallow minimum at 
ksT ^ A (Fig. The behaviour in the empty orbital regime (eo/A = +1) is similar to 



that in the mixed valent regime except that the shallow minimum at ksT eo ~ A becomes 
a shoulder (Fig. [171) . 

This complicated behaviour of the thermopower has a clear interpretation in terms of 
the temperature dependence of the spectral densities. We consider first the Kondo regime. 
The low temperature behaviour T « of the thermopower obtained from the Sommer- 
feld expansion S{T) ■yT shows that S is positive and large due to the strongly enhanced 
values of 7 in the Kondo regime (7 ~ 1/Tk)- At temperatures T of order Tk for which the 
thermopower is determined by the Kondo resonance, we see from (|^) that the sign of the 
thermopower depends on the slope of the spectral density at the Fermi level ( ^poCi^'^) ] 
This slope decreases and eventually changes sign for T > Tk as a result of the strong tem- 
perature dependence of the Kondo resonance. This sign change at T ^ Tk is characteristic 
of the Kondo regime. As the temperature is increased to ksT ^ A, charge fluctuations be- 
come increasingly more important and there is a transfer of spectral weight from the lower 
satellite peak at eo to the upper one at eo + U (Fig. |^). Holes are scattered increasingly more 
than electrons as the temperature is increased and therefore the thermopower increases for 
ksT > A. Hence a broad minimum arises at fc^T ^ A. 

The behaviour of the thermopower in the mixed valent and empty orbital regimes can be 
explained in a similar way in terms of the temperature dependence of the spectral density. In 
the mixed valent regime the interactions renormalize the bare local level eo to lie at e~o above 
the Fermi level (see Fig. ^). The thermopower can be analyzed in terms of the resonance 
at eo of width A in a similar way to the above analysis for the Kondo resonance. The 
low temperature maximum occurs at ksT ^ A/3 but the low temperature enhancement 
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S{T) ~ 7T is only of order 7 ~ 1/A (neglecting the phase factor cot(7mo/2)) . There is 
a maximum at ksT ^ A, beyond which the thermopower decreases. In this case however, 
the resonance above the Fermi level does not disappear but is increasingly broadened with 
temperature (see Fig. Hence the thermopower does not change sign. In the empty 
orbital regime, the behaviour of the thermopower can be similarly explained in terms of the 
temperature dependence of the spectral density. 

We now consider the effect of non-resonant scattering of conduction electrons on the 
thermopower. The effects of non-resonant scattering are most dramatic for the symmetric 
case, where in the absence of such scattering the thermopower vanishes identically. The 
presence of non-resonant scattering leads to interference effects between the resonant / = 
channel and non-resonant / = 1 channels making the thermopower finite even for the 
symmetric case (see Appendix 1 for details of how these interference terms arise). This is 
seen explicitly at low temperatures by applying the Sommerfeld expansion to (^ keeping 
the full transport time with inclusion of non-resonant terms, 

sin(2?7o(ei;^) - 2r]i) 



\e\ [sm (?7o(eF) - Vi) + Pr. 
Here, tiq^u) is the resonant phase shift. 



+ 0{T^). (46) 



and rji is the phase-shift of conduction electrons scattering in the / = 1 channel (p„ contains 
the effects of higher angular momentum scattering and was discussed earlier). For rjolep) = 
7r/2 corresponding to the symmetric case the above expression gives a finite thermopower 
due to the non-resonant terms. 



In Fig. |T8| the thermopower for the symmetric case and several values of the non-resonant 
phase shift rji is shown. The thermopower can be positive or negative at low temperature 
depending on whether rji is negative or positive. This is clear since a positive r/i corresponds 
to a repulsive potential which decreases the charge on the impurity making the thermopower 
the same as for less than half-filling, i.e. positive at low temperature and conversely for t]i 
negative. 
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The effect of non-resonant scattering of conduction electrons on the thermopower in the 



asymmetric cases is shown in Fig. |T9] and Fig. for the Kondo and mixed valent cases 
respectively. In the mixed valent case we see that a small positive non-resonant phase shift 
can change the minimum at ksT ^ A (Fig. ^) into a shoulder making the thermopower 
resemble that of the empty orbital case. On the other hand a small negative non-resonant 
phase shift accentuates the minimum aX ksT /\, and makes the thermopower more similar 
to that of the Kondo case. It is possible that a sufficiently large negative phase shift rji 
could make the thermopower in the mixed valent case change sign, but in this case, which 
corresponds to strong non-resonant scattering, the frequency dependence of r/i should be 
taken into account. In the Kondo regime similar trends are found on varying the non- 
resonant phase shift from negative to positive values. In the weak Kondo regime, cq/A = —2 
(Fig. |T9|), for which = 0.18 A, a sufficiently large value for rji > 0.042 can change the 
sign of the high energy minimum at ^ A making the thermopower resemble that of 
the mixed valent case. The effects of non-resonant scattering can be understood as for 
the symmetric case as arising from a change in the local screening charge on the impurity. 
A difference in the thermopower of Fig. jTBI for the symmetric case and Fig. |TU| for the 
asymmetric case is that the extremum at ksT ^ A in the symmetric case is very much 
suppressed. This is expected since this extremum is associated with the charge fluctuations 
which in the symmetric case are strongly suppressed and only included in the present case as 
a result of the non-resonant scattering terms. We see then that the thermopower is sensitive 
to small changes in the local environment of the magnetic ion which has been modelled here 
by including non-resonant scattering of conduction electrons. The sensitivity arises because 
of interference between the / = 1 channel and the resonant / = channel, in contrast to 
the I > 1 channels which have negligible influence since they do not couple to the resonant 
channel. 
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Thermal Conductivity and Lorenz number 



The thermal conductivity is shown in Fig. 21. It shows the expected Fermi hquid be- 
haviour K,{T)/aT ~ 1 + I3{T /TxY low temperature (T << Tk) with hnear coefficient 
a ~ T^tp) ~ 1/ sin^(7mo/2) as shown in Fig. 22 and in accordance with the Friedel sum 
rule. The coefficient /3, i.e. the coefficient of k(T), has also been extracted from 
the numerical results by using a least squares fit in the region < T /Tk < 0.01. The 
behaviour is shown in Fig. ^ and the results extracted for (3 are listed in Table |. It lies 
within 15-25 % of the exact result /? = 15.83 calculated above using Fermi liquid theory. 
The agreement is surprisingly good since we are dealing with corrections to k(T) and the 
extraction of (3 also relies on an accurate extraction of the linear coefficient a. The resis- 
tivity coefficient comes from a term and is consequently more accurate. A distinct low 
temperature anomaly is absent from /t(T), and its behaviour is similar in both the Kondo 
and mixed valent regimes. Universality in k(T)/T extends up to at least lOT/^ in the Kondo 
regime. The Lorenz number ratio L/Lq where L = np/T and Lq = vr^/cs/Se is shown in 
Fig. |2^. Deviations from the Wiedemann-Franz law are found for T > O.IT;^. In the Kondo 
regime L{T), like k,{T), is universal up to approximately lOT^;. In the Kondo regime it also 
exhibits a maximum at T ^ whereas in the mixed valent there is only a shoulder at 
^ A, and similarly in the empty orbital regime. The inclusion of weak non-resonant 
scattering has little effect on k(T). 

Hall coefficient 



Fig. ^ shows the temperature dependence of the Hall coefficient. This quantity exhibits 
a low temperature peak at T Tk in the Kondo regime but not in the mixed valent or 
empty orbital regimes. Universality for this quantity persists only up to T Tk- Here we 
have only considered the ordinary Hall effect and have shown that this exhibits an anomaly 
at T Tk- A more important contribution to the Hall effect arises from skew scattering 
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of conduction electrons. Up spin electrons are scattered differently in a magnetic field than 
down spin electrons as a result of the splitting of the Kondo resonance by the magnetic field. 
This contribution should be much more important than that coming from the Lorenz term 
considered here, but its treatment is outside the scope of the present work. 



DISCUSSION AND CONCLUSIONS 

In this paper we have given a detailed description of the generalization of the numeri- 
cal renormalization group to the calculation of finite temperature Green's functions of the 
Anderson model. The method was applied to obtain the single-particle spectral densities 
and transport properties of the Anderson model in the strong correlation limit fZ/vrA >> 1 
for a range of temperatures and parameter regimes, as well as the occupation number and 
specific heat. We distinguish three parameter regimes in the limit of strong correlation. The 
Kondo regime, eo << —A, uq ^ 1, characterized by a low energy scale ksTo = ksTx, the 
valence fluctuation regime, —A < cq < 0, 0.3 < r^o ^ 0.8, characterized by a low energy 
scale fcsTo = A, and the empty orbital regime, eo > A, uq < 0.3, characterized by an energy 
scale ksTQ ~ cq ~ cq. We also distinguish three temperature ranges, T « Tq corresponding 
to the Fermi liquid regime, T » Tq corresponding to the high temperature perturbative 
regime and T ^ Tq corresponding to the cross-over region. A consistent picture of the 
thermodynamics, spectral densities and transport coefficients emerges from the numerical 
renormalization group calculations in these different regimes. 

The calculations of the single particle excitation spectrum show that at low temperatures 
the spectral function in the Kondo regime is characterized by a Kondo resonance of negligible 
spectral weight, centered around ep, and two atomic like resonances at eo and eo + U, which 
carry most of the spectral weight. The one below the Fermi level at eo contains most of the 
screening charge while the one above at eo + t/ is empty. This picture is well known and 



emerged very clearly in the finite order perturbation calculations of Yamada |31|. Using the 



numerical renormalization technique we have been able to calculate it more accurately and 
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obtain the correct exponential dependence of the width Tk on eo and U and also describe 
the high temperature regime where there are corrections that depend logarithmically on T. 

The spectral function for valence fiuctuators is characterized by two peaks only: a par- 
tially filled one at the renormalized position eg > ep and an empty one at eio + U . In the 
empty orbital regime the low temperature spectral function in the strong coupling limit 
is just an empty virtual bound state. The characteristic low energy scale for the spectral 
function is set by Tq. These features of the spectral density are fully consistent with the 
Fermi liquid groundstate as revealed by the thermodyamics [§] . The Friedel sum rule which 
is a good test on the accuracy of the method at both high and low energies was found to be 
satisfied to within 3% accuracy in all parameter regimes. Increase of temperature leads to 
a renormalization of the excitation spectrum and for T > Tq the spectral function changes 
in a qualitative way 0. At high temperatures the spectral function is characterized by two 
atomic like peaks separated in energy by U, consistent with the thermodynamic properties 
which exhibit local moment behaviour. In the mixed valent and empty orbital regimes the 
higher peak is strongly suppressed. 

The transport coefficients for the Anderson model can be related to the transport re- 
laxation time by using linear response theory. This can be further related to the single 
particle spectral function which allows us to discuss the transport coefficients in terms of 
elementary excitations in the resonant channel. The characteristic low energy scale of the 
spectral function, Tq, defines the temperature scale for transport properties. 

In the Kondo regime at low temperature, T << Tq = T^, we showed that the trans- 
port coefficients agreed with the predictions of Fermi liquid theory. For T « T^, the 
localized spins are screened and the transport currents relax via excitations in the local 
Fermi liquid state formed by the Kondo screening mechanism. The transport coefficients 
were characterized by power law dependencies on T /Tk (|3^-p^ with universal coefficients 
for appropriately scaled quantities. The numerical results for the term of the resistivity 
and the T^ term of the thermal conductivity in the symmetric case agreed remarkably well 
with the analytic values deduced from Fermi liquid theory. With increasing asymmetry the 
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numerical results indicate that these coefficients change but only slightly for no > 0.8. The 
thermopower shows universal behaviour for T < when rescaled by its zero temperature 
slope Universality was found to extend up to at least IOTr- for the resistivity and 

thermal conductivity but only up to Tk for the Hall coefficient. The transport properties in 
this regime are understood in terms of excitations in the strongly renormalized Fermi liquid 
ground state. On increasing temperature through the cross-over region T ^ Tk logarithmic 
terms become important for the resistivity and these are clearly seen for T > Tk- The 
thermopower changes sign, a result of the disappearance for T > Tk of the Kondo reso- 
nance in the single-particle spectral density. There is also a minimum in the Hall coefficient 
at T Tk- In the high temperature perturbative regime, T >> Tk, the spectral density 
consists of two peaks, the occupancy uq approaches 1 and the transport properties are those 
of electrons coupled weakly to local spins. The thermopower develops a broad minimum at 
fc^T ^ A due to a shift of spectral weight from the peak below the Fermi level to that above 
as charge fluctuations become important. 

In the mixed valent and empty orbital regimes the transport coefficients at low tem- 
perature T << To again reflect the relaxation of the transport currents via the excitations 
in the renormalized local Fermi liquid ground state. However, since A >> ksTK, these 
renormalizations of order I/Tq are smaller than for the Kondo case. The low temperature 
transport coefficients again have power law dependencies on T/Tq, but the coefficients are 
non-universal and depend strongly on the occupation of the localized level no. 

Finally, we investigated the effects of non-resonant scattering of conduction electrons 
on the temperature dependence of the transport properties. A general expression for the 
transport time in the presence of such scattering was derived which shows that the result- 
ing interference effects between the resonant and non-resonant channels are not negligible. 
Hence, to relate the transport coefficients of the Anderson model to the experimental data 
on dilute alloys one should take into account the non-resonant scattering. For example, the 
overall shape of the thermopower curves calculated for the Anderson model with two scat- 
tering channels resembles quite strikingly the experimental data. The interference effects 

28 



between resonant and non-resonant scattering channels led to a strong dependence of the 
sign of the thermopower on the non-resonant phase shift. This type of effect could be im- 
portant for systems with no close to 1 or for systems with a Kondo resonance approximately 
symmetric about the Fermi level. Doping such systems with non-magnetic impurities could 
lead to the above type of interference effects making the sign of the thermopower sensitive 



to the local environment. Such effects have been observed In the symmetric Kondo 
regime the inclusion of non-resonant scattering was important to obtain the non-vanishing 
enhanced thermopowers for half-filled systems (the no = 1). 

To summarize, we have shown that the numerical renormalization group method allows a 
unified description of the transport, thermodynamic and spectral properties of the Anderson 
model in the strong correlation limit. The different features below and above the character- 
istic scale To reflect differences in the structure of the excitation spectrum at high and low 
temperatures. The NRG method presented here could be used to provide further insights 
into the behaviour of the Anderson model. For instance, the results for the spectral density 
could be used to deduce the spin and charge densities induced in the conduction electrons at 
arbitrary distances from the impurity |^ . The interference effects for the transport proper- 
ties could also be studied for two Kondo (spin 1/2) impurities or for the case of the Kondo 
effect in the / = channel and a virtual bound state in the / = 1 channel. This would be 
relevant to rare earth systems in which there is an excited crystal electric field level close 
to the Kondo ground state. There are many other similarly interesting physical situations 
which occur in particular materials. They could be modelled and a broad range of physical 
properties calculated via this one technique to compare with experiment. 
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APPENDIX A: DERIVATION OF THE TRANSPORT TIME IN THE PRESENCE 

OF NON RESONANT SCATTERING 

The transport time for scattering from a dilute concentration, Ui << 1, of magnetic im- 
purities described by the Hamiltonian (|1|) can be obtained by examining the Kubo expression 
for the conductivity o"(T), 

air) = -lim^^o- Im Uretitu), (Al) 

UJ 

where the current-current correlation function Il{iuj) = —-^ Jo e^'^^ < TT-j(r)j(0) > and 
the current operator is that for free-electrons j = — — Z^ko-k cj^o-'^ko-- The diagrammatic 



representation of Il{iuJn) is shown in Fig. 26. In the absence of non-resonant scattering the 



vertex corrections vanish identically [jl3| in the dilute limit leaving just the (dressed) bubble 



diagram Ho. In terms of the dressed conduction electron Green's function G(k, u;) this is 
noM = ^^J2^'WGik,iuJn + tto)Gik,tUn). (A2) 

Here we evaluate the conductivity including the vertex corrections since these give rise to 
interference terms between the / = and / = 1 channels. The conductivity is given by 

a{T) = ^ ( ^ ) [P(e -i5,e + i5) - Re[P{e + i5,e + i5)]] de, (A3) 
In J-oo \ oe J 

where the singular part behaving like l/n^ comes from the first term with P(e + i6,e — i6) 



28| . The latter is given in terms of the conduction electron Green's function G(k,e) and 



vector vertex function F by 

2e2 1 

P(e -iS,e + iS) = --y G(k, e + i5)G(k, e - i5)k.T(k, e-i5,e + i5). (A4) 

3m^ V ^ 

The vertex function T satisfies the following equation 
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r(k, e - e + = k + / — — r(k', e-i5,e + i5)W^^i{e -i5,e + ^5)G(k', e - i5)G{e + i5) 

(A5) 

where W is the irreducible vertex. 

To solve this equation we introduce a scalar vertex function 71 defined by r(k, e — e + 
i5) = k 7i(k, e) and note that the Green's function product in ( |X5| ) is just A(k, e)/|S'^(k, e)| 
where A(k,e) is the spectral density and E^(k, e) the imaginary part of the self-energy of 
the conduction electrons. In the dilute limit, Ui 0, we also have that W^kk'(^ — i6,e + 



i5) — > ?T,jTkk'(e — «5)Tk/k(e + i5), where T is the energy dependent T-matrix, |S'^(k', ek/)| — »• 
nj|JmTk'k'| and v4(k', e) 27r5(e — Ck')- The resulting equation for the scalar vertex function 
after some algebra becomes 

/IN T , ^ f ^ M' /o sin d' cose'6{e - ek')Tkk'(e - iS)T]^>]^{e + i6) ji{k', e) 

7i k,e = 1 + — / dek'de V2mek' — . A6 

Air J |imik'k'(^ + 

where 6' defined by cos(^^') = k.k'/k^ is the angle between incoming and scattered k states 
and a free electron density of states has been used to rewrite the integrals. The T-matrix 
in the denominator of ( [A6[ ) is for forward scattering only so it is independent of 6' and can 
therefore be taken outside the integral and evaluated at k' = k (note the 6 function in ( |A6D ) . 
The scalar vertex function 71 (k', e) is also independent of angle and can also be taken outside 
the integral (evaluated at k' = k). Hence the closed expression, 

7i(fc, e) = l + "^TT^f^^—^MK e), (AT) 
47r \ImTkk[ek + m)| 



is obtained with solution, 



N \ImTkk{ek + , . 

\ImTkk{ek + to)\ - -j-h 



In the above expressions Ik given by 

4 = r de' sin e' cos ^'Tkk' (ck - «5)Tk'k(ek + i^) (A9) 
Jo 



contains the interference effects. 
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Substituting the scalar vertex function 71 from ( |X8| ) into the expression for P{e—iS, e+iS) 
(using the previosuly quoted expression for the Green's function product) gives 

P(e _ ,5, e + ^5) = ^ fcM(k, e)7i(fc, e)/S^(k, e) (AlO) 

k 



and the latter into the conductivity formula (|A1|) gives for the transport time r the expression 
T{ek,T)-' = 2ni{\Imnk{ek + -"^h). (All) 

The first part of this expression is just the usual result for s-wave scattering only and the 
second part is the contribution coming from interference effects between the resonant / = 
channel and the non-resonant channels. To evaluate (^^) and ( [All| ) we need the T-matrix. 



This can be obtained for the Hamiltonian (|I|) by aplying the equation of motion method to 
the conduction electron Green's function CkkV =<< Cko-; c^/^. >>• This gives 

G'kkv(^ + ^S) = G°kv(^ + + G^kvl^ + ^5)Tkk'(^ + ^S)GL'A^ + (A12) 

where 

^kk'(^^ + iS) = V"okG'o(a; + iS)Vok> + ^ VimkGi{uj + iS)Vimk', (A13) 

l=lm 

^kk'o- unperturbed conduction electron Green's function and Gi,l = 0, 1, . . . are the 

local Green's functions for the Z = 0, 1, . . . channels respectively. The hybridization matrix 
elements VJ^k have the angular dependence VJ^k =< k\V\l,m >~ VikYimi^k) so the T 
matrix has the form: 

TkkKc^ + ^5) = ^(2/ + l)Pi{cosie'))VikGiiuj + t6)Vik^ (A14) 
1=0 

= VokGoiio + t6)Vok' + 3cos{e')VikGi{uj + t6)Vik' + • • • (A15) 

(A16) 

where EmYiU^k)Yi^{nk') ~ (2/ + 1) Pi{cos{9')) has been used and 9' = k.k'/k^. The T 
matrix has an explicit dependence on the angle 6' which gives rise to interference effects 
between the resonant and non-resonant channels when substituted into the expression for 
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the transport time ( |A11|) . These interference effects arise in the expression for Ik as cross- 
terms involving GqGi + GqG\. On carrying out the 9' integrations in ( |A9| ) and expressing 
the non-resonant Green's function Gi{uj) in terms of a corresponding phase shift tii{ijj) we 
arrive at the formula for the transport time of the multi-channel Anderson impurity model: 

1 1 



Re Go{u;,T) . 

cos 2r]i - — — — sm 2r/i 

Im Go{uj, I ) 



+ Pn. (A17) 



T(a;,T) T%uj,T) 

The resonant part of the transport time, ro(c<j, T) is given by the usual T-matrix expression 
in the absence of the interference effects, 

— ^ = Apo(^, T) = -A Im Goioo, T)/7r, (A18) 

where, A is the unrenormalized width of the / = resonance and Gq{uj,T) is the / = 
Green's function, which describes the many-body effects. The analysis here was for the 
electrical conductivity. The same transport time ( |A17|) is found by repeating this analysis 
for the other current-current correlation functions (In the above expressions for the inverse 
transport times, a constant factor 2ni/N{0), where A^(0) is the unperturbed conduction 
electron density of states at the Fermi level, has been omitted) . 



APPENDIX B: THE NUMERICAL RENORMALIZATION GROUP METHOD 

The following continuum version for the resonant part of (0) is useful for introducing the 
logarithmic discretization approximation 

H = eocl^coa + Uno^noi + Vq J ^ dk{cl^Coa + cl^Cka) + J ^ kc{^Ckadk. (Bl) 

In obtaining this from the resonant part of (|l]), a partial wave expansion for the con- 
duction states around the impurity has been used has been used and the Hamiltonian 
has been written in the energy representation. The approximation has been made that 
A(cij) = TiY^k |Vbfcp5(u; — Cfe) = A is independent of frequency so that only s-wave partial 
wave states couple to the impurity (with strength Vq = J A/'KN{eF))- AH energies are 
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measured relative to the half bandwidth D = 1 and we use the convention that repeated 
spin indices are summed over, (full details of the above may be found in 0). The contin- 
uum conduction band [—1, 1] is now discretized as shown in Fig. and a new set of basis 
states for the conduction electrons is introduced in each interval ±[A~^"+^\ A^"] using the 
following wavefunctions 

f nn^Mi/a e^'^"^'' for A"""! <±k< A-" 
V^i(^) = ""^'-'/^^ ' (B2) 
otherwise 

Here p is a Fourier harmonic index and Un = 27rA"'/(l — A^^). The operators Cka can then 
be expanded in terms of new operators Qnpa, bnpa labelled by the interval n and the harmonic 
index p 



Cka = [(^npai'tp + bnpai'np] ■ (B3) 
np 

In terms of these operators, the Hamiltonian ( pi|) becomes 0|, 
H = eocl^coa + Uno^noi 

oo 
n=0 

^ np 

+ E E A-"(at,,a„,v - bl^b^p,^)ef^ . (B4) 

n p=^p' 

The coupling of the local state Cqo- to just the p = harmonic is a consequence of the 
assumption that the hybridization matrix elements are independent of energy. Hence the 
conduction electron orbitals a^p, bnp for p 7^ only couple to the impurity state cqo- indirectly 
via their coupling to the anQ,bno in the last term of ( |B4| ) . This coupling is weak, being 
proportional to (1 — A^^), and vanishes in the continuum limit A — > 1, so these states 
may be expected to contribute little to the impurity properties compared to the p = 
states. This is indeed the case as can be seen in Fig. |2^, where the probability density 
in real space of the wave-packet states ipnpik) is shown. It can be seen that the p > 
orbitals are localized at increasing distances from the impurity, whereas the p = orbitals 
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overlap strongly with the impurity. This together with their weak coupling to the impurity 
ensures that they can be neglected in calculating local impurity properties. This logarithmic 
discretization approximation has been shown to give rapidly convergent results (as A — > 1) 
for thermodynamic averages. The deviations away from the continuum limit, A = 1, being 
proportional to e '"(^) ^'^^ approximation amounts to approximating the continuum 

conduction band Hamiltonian He, by a discrete one: 

^ hi + A-i) ^ A-"(4o.a„oa - blo^bnoa). (B5) 

^ n 

Consequently a discrete approximation to the Anderson model is obtained 

OD 

n=0 

+ l(l + A-i)EA-"(aLa„.-6iA.), (B6) 

^ n 

where the p index has now been dropped since only p = states are included. This discrete 
Hamiltonian is more convenient for a numerical treatment. 

The next step is to convert the discrete Hamiltonian ( P6|) into a linear chain form suitable 
for an iterative diagonalization starting from a small chain and diagonalizing sucessively 
longer chains by adding a site or energy shell at each stage. If we define a new local 
conduction electron orbital by |Ocr >= /Q^|vac > where 

-| oo 

fl = [^(1 + E A-"/^((aL + O, (B7) 

^ n=0 

then the Lanczos algorithm can be used with starting vector |Ocr > to tri-diagonalize 

oo oo 

-f^c = E ^nflafna + E ^nifn+lafna + fnafn+la)) (B8) 
n=0 n=0 

where en,\n,n = 0,1,... are the site and hopping energies respectively. For the half- 
filled Fermi sea |vac > the e„ are all zero and Wilson has obtained the expression A„ = 
|(1 + A^^)A^"/^^„, where ^„ = ^n(A) 1 for n >> 1. This converts ( [B6| ) into a semi- 
infinite chain with the impurity at the origin coupled to the local conduction state \0a > 
via the hybridization Vq (re-scaled to keep |0o" > normalized) 
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H = eocl^coa + f/rioT^oi + K)(/oVcoa + H.c.) 

1 oo 

+ -(1 + Y: UfUlJn. + flfn+l.) (B9) 

^ n=0 

i/^ = A(^-l)/2[ifo + H,y, + Y: A-"/2en(/I+l./„. + /L/n+1.)]. (BIO) 

n=0 

In ( piO|) we have defined a sequence of rescaled Hamiltonians Hn whose smallest term is 
always of order 1 (the bare parameters eo, U, Vq are also rescaled by the factor |(1 + A^^)). 
The procedure is now to diagonalize this sequence of Hamiltonians iteratively and extract 
the spectrum and eigenstates on successively lower energy scales cjat ~ A-(^-i)/2. The 
procedure starts by diagonalizing the impurity part, Hq = coCqo-Coo- + UriQ-^nQi of (|E^ ) and 
then adding the coupling to the local conduction electron orbital Vo(/oo-Cocr + H.c). Suc- 
cessive energy shells ^nifl+icrfna + flafn+ia), n = 0,1, . . . are then added and the resulting 
Hamiltonian diagonalized to give the many-body eigenvalues and eigenvectors |p >Ar 
of the corresponding shell or cluster. The total number of electrons, A'^g, total spin, S, 
and z-component of total spin, Sz, are conserved quantities and can be used to label the 
eigenstates together with an index r = 1, R^^s where R^^s dimension of the subspace 

(NeS). If Hn has been diagonalized, Hn = Y.N,ss.r EN,s^N,ss,r,NeSS,r, then from ( P10|) 
the matrix for H^^i can be obtained from 

The basis used for Hn+i is the product basis > of eigenstates, \p >= {NeSS^r >, of 
Hn, and states \i > from site + 1 (i.e. \i >= |0 >, | t>) I i>) I Ti>)- In this basis the 
matrix elements of i^Ar+i are 

<p,i\HN+l\q,j > = A^^^Sp^gSijE^ + ^N{<P,i\fN+lafNa\qJ > + < P,i\fNafN+la\qJ >) 

+ e7v(-l)^'=-(< i\fh+lJj X p\fNa\q > + < PifLk X AfN+la\j >), (B12) 

where fN+ia have been commuted past the eigenstates \p >, \q > with a sign change (— 1)^=''' 
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depending on the number of electrons, N^^q = Ne^p, in these states. Since is already 
diagonalized the matrix elements < p\fNa\(l > can be calculated as follows. The unitary 
transformation between the eigenstates of if at and the product basis is known and is given 
by 

\p>N=J2UN{p,ri)\r>N-i\i>, (B13) 

r,i 

where UN{p,ri) is the matrix of eigenvectors of H^, \r >n-i denotes an eigenvector of the 
previous cluster (A^ — 1) and |i > is one of the four states given above, but for site (not 
N + 1). Hence the required matrix elements are given by 

< P\fNa\q >N =J2J2u;f(P^ri)UN{p,r'i') < i\ < r\fNa\r' > \i' > 

r,i r',i' 

= T.T.UUP,r^)UN{q,r'^')6ry{-lf^- < ^\fNa\^' >, (B14) 



which involves only known expressions. The matrix elements < i\fN+ia\j > in ( |B12| ) are also 



known so the matrix of H^^i can be set up in terms of the eigenvalues and matrix elements 
in ( |B14| ). Hjy+i can then be diagonalized and the procedure repeated to obtain the spectrum 

iV— 1 

on successively lower energy scales, ujn^2,^n+3, ■ ■ where ujn ~ -DA 2- is the smallest 
scale in Hj^. In practice since the number the number of many-body states in Hi<^ grows like 
4^ it is not possible to retain all states after about N = 7 (even after symmetry has been 
used to reduce the size of the matrices). For N > 7 only the lowest 1000 or so states of Hjsi 
are retained. The truncation of the spectrum in Hj^r restricts the reliable range of excitations 
uj to uj]\f < uj < Kuj]\f where K is a constant dependent on A {K ^ 10 for A = 3). The lower 
excitations are calculated more reliably in successive iterations, whilst information on the 
higher excitations is contained in previous iterations. The eigenvalues are used to calculate 
the partition function Zjy^Pjy = l/fc^T/v) and free energy FiTjy) = — fc^ ln(^Af(^iv)) for the 
total system at a decreasing sequence of temperatures T/v = uj^/kB ~ DhT^^^^'^/'^ ,N = 
0, 1, . . .. The thermodynamic properties are then extracted by first subtracting out the 
conduction electron contribution to the free energy, Fc(T/v) The specific heat is then 

given by the second derivative of the impurity free energy, Fimp(T) = F(T) — Fc(T), 
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(B15) 
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FIG. 1. The local level ocuupancy no{T) over the whole temperature range in various regimes. 

FIG. 2. The specific heat in units of over the whole temperature range and in the Kondo, 
mixed valent and empty orbital regimes. 

FIG. 3. The specific heat exhibits Fermi liquid behaviour, C{T) = ^T, at low temperature, 
T « Tk, in all parameter regimes. Note the very much enhanced values of 7 in the Kondo cases 
compared to the mixed valent and empty orbital case. 

FIG. 4. The specific heat in the Kondo regime, showing the universal behaviour for T < Tk- 
However, when the local level position is well separated from the Fermi level, universality extends up 
to almost WTk over almost the entire range of the low energy peak. The two disticnt energy scales, 
Tk, governing the low energy peak and A >> Tk governing the high energy charge fluctuation 
peak are clearly evident. 

FIG. 5. The T = spectral density in the Kondo regime for several positions of the local level: 
the lower satellite peak is at the local level position eo. This takes values cq = —U/2 (symmetric 
case), eo/A = —4, —3, —2 (asymmetric cases). For eo/A = —2 the lower satellite peak forms a 
shoulder and has partly merged with the Kondo resonance. The high energy satellite peak is 
separated from the lower one by the Coulomb energy U. The inset shows the Kondo resonance in 
more detail. The Kondo resonance becomes broader as — eo decreases. 

FIG. 6. The finite temperature spectral density for the symmetric case. The inset is for p(u}, T) 
plotted versus T/Tk and shows in more detail the region around the Fermi level. The satellite 
peaks acquire some temperature dependence for T A. 

FIG. 7. The spectral densities for the asymmetric case in the Kondo regime. At high temper- 
atures T > A there is a shift of spectral weight from the lower to the upper satellite peak. The 
Kondo resonance, shown more clearly in the inset where it is plotted versus T/Tk, disappears for 
T » Tk. 
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FIG. 8. The T = spectral densities in the mixed valent (eo/A = —1,0) and empty orbital 
(eo/A = +1) regimes. The resonant level of width of order A is renormalized by the interactions 
to lie above the Fermi level. There is no Kondo resonance in these cases and the upper satellite 
peak at approximately eo + U has very little weight. 

FIG. 9. The temperature dependence of the spectral densities in the mixed valent regime. The 
renormalized resonant level of width A lying at eo > acquires a strong temperature dependence 

on a scale of T ~ A. 

FIG. 10. The temperature dependence of the spectral densities in the empty orbital regime. 
The relevant scale is the renormalized level eo > which for the present choice of parameters is 
approximately A. On a scale of eo ~ A this resonance has a strong dependence on temperature. 

FIG. 11. The electrical resistivity in various parameter regimes and over the whole temperature 
range. The curves are plotted versus the reduced temperature T/Tq, where To = A in the mixed 

valent and empty orbital regimes and Tq = Tk in the Kondo regime. The zero temperature value, 
p{T = 0), is fixed by the Friedel sum rule to be sin2(7rno(T = 0)/2)/7rA. 

FIG. 12. The resistivity in the symmetric case with a fit of the Hamann expression with 
Tkh = Tk/1-2. The Hamann result fits the high temperature range, and is very close to the NRG 
data for Tk <T< WTk- It fails at low and very high temperatures. 

FIG. 13. The scaled resistivity in the Kondo regime showing the universal behaviour at low 
temperature up to approximately 5Tk- The inset for 1 — (p(r)/p(0)) versus {T/TkY shows the 
expected Fermi liquid behaviour for the resistivity at low temperature T < O.ITk- The coefficient, 
c, of the term in the resistivity is found to lie within 8 % of the exact result in all cases. 

FIG. 14. The influence of non-resonant scattering of conduction electrons on the resistivity in 
the symmetric Kondo regime. 
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FIG. 15. S{T)/{T/Tk) showing the hnear in T Fermi hquid behaviour of the thermopower 
for T « Tk- The different curves are for local level positions eo = — 4A(o), eo = — 3A(n), 
eo = — 2A(o), eo = — 1A(A), eo = 0(v)5 eo = The inset shows the first three cases in 

more detail. 

FIG. 16. The thermopower iS'(r) in the Kondo regime. The low energy maximum is at a 
temperature T in the range T^/S <T< Tk whilst the broad minimum at higher temperature is 
at T Ri A. 

FIG. 17. The thermopower S{T) in the mixed valent and empty orbital regimes. The low 
energy maximum in the mixed valent case is in the range A/3 < T < A and the minimum at 
higher temperature is at T A. In the empty orbital case, there is only a shoulder at T A 

FIG. 18. The thermopower >S'(r) in the symmetric case with inclusion of non-resonant scat- 
tering. The low energy minimum is at T Ri Tk and the extremum at higher temperature is at 
T Ri 50Tk Ri 0.9A. 

FIG. 19. The effect of including non-resonant scattering on the thermopower in the Kondo 
regime eo = — 2A. For a large enough non-resonant phase shift rji > 0.0427r the broad minimum 
at r Ri A changes sign. The low energy maximum is at T Ri Tk/S. 

FIG. 20. The effect of including non-resonant scattering on the thermopower <S'(r) in the mixed 
valent regime eo = 0. There is only a broad maximum at T Ri A/3 which disappears when the 
non-resonant phase-shift rji is increased. 

FIG. 21. The thermal conductivity l^l^lj'yT)o Plotted versus 7 T. The solid lines with symbols are 
for the Kondo regime. The mixed valent cases are eo = — A (long dashed) and eo = (dashed), and 
eo = —A (dotted) is the empty orbital case. (The single point for the symmetric case which falls off 
the universal curve is due insufficient accuracy of the integrations at that particular temperature) 
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FIG. 22. The linear coefficient of the thermal conductivity a = (k(T)/T)o versus ^.^^j^^^^^^^ 
(see Table I for the values). 

FIG. 23. The coefficient of K{T)/a T in the Kondo regime (see Table I for the values). 

FIG. 24. The Lorenz number ratio L(T)/Lq where Lq = Tr^fes/Se versus 'jT. 

FIG. 25. The Hall coefficient R{T)/Rinf versus 7 T where i?j„/ = —l/n\e\c. 

FIG. 26. The diagrammatic representation of the current— current correlation function Il{iu>). 
The solid lines represent full conduction electron Green's functions and F is the two-particel scat- 
tering vertex. 

FIG. 27. The conduction band [— 1,+1] logarithmically discretized with discretization param- 
eter A > 1 

FIG. 28. The normalized conduction electron orbitals Ikpripupif)]'^ as a function of distance 
kpr from the impurity for n = 5 and p = 0,1, ... ,5. 
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TABLE I. The first column shows the local level positions used in the calculations, and the other 
columns show various data extracted, p^"^^'^'' = SE-^^aM^ where no is the local level occupancy 
at T = and is extracted from the partition function. The coefficient of the resistivity in 
the Kondo regime is defined by c = limT-^-o ( ^~(T/Ti)'^^^ ) ' which according to Fermi liquid theory 
should be 6.088 in the symmetric case with small corrections in the asymmetric Kondo regime. 



The linear coefficient of the thermal conductivity is a = ("T^)^ ^ P is the coefficient of 
K,{T)/aT which should be 15.83 in the symmetric case. 





no 


7/fe| 


Tk/A 


Po 




% error 


c 


% error 


a 


/3 


.o = -f 


1.00 


8055.5 


0.0180 


31.83 


32.61 


-2.4% 


5.7 


-6.4% 


0.102 


11.6 


eo = -4A 


0.93 


4224.0 


0.0.346 


31.45 


32.13 


-2.1% 


5.8 


-4.7% 


0.101 


13.5 


eo = -3A 


0.88 


2069.0 


0.0690 


31.26 


31.33 


-0.2% 


6.4 


+5.1% 


0.106 


12.6 


eo = -2A 


0.78 


832.0 


0.1790 


28.34 


28.93 


-2.0% 


6.6 


+8A% 


0.113 


13.2 


eo = -A 


0.63 


281.0 


0.5900 


22.14 


22.32 


-0.8% 






0.147 




eo =0 


0.45 


100.7 


2.5100 


13.18 


13.36 


-1.4% 






0.246 




eo = +A 


0.31 


42.2 


13.660 


6.85 


6.79 


-1-0.9% 






0.496 
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